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Timing is crucial to many aspects of human performance. To better understand its neural 
underpinnings, we used event-related fMRI to examine the time course of activation associated with 
different components of a time perception task. We distinguished systems associated with encoding 
time intervals from those related to comparing intervals and implementing a response. Activation in 
the basal ganglia occurred early, and was uniquely associated with encoding time intervals, whereas 
cerebellar activation unfolded late, suggesting an involvement in processes other than explicit 
timing. Early cortical activation associated with encoding of time intervals was observed in the right 
inferior parietal cortex and bilateral premotor cortex, implicating these systems in attention and 
temporary maintenance of intervals. Late activation in the right dorsolateral prefrontal cortex 
emerged during comparison of time intervals. Our results illustrate a dynamic network of cortical- 
subcortical activation associated with different components of temporal information processing. 


Humans are remarkably proficient at perceiving the passage of 
time and producing precisely timed behaviors, many of which 
depend upon explicit prospective temporal judgments. For these 
events, multiple processes seem to determine our subjective per- 
ception of current time for intervals lasting several hundreds of 
milliseconds to several seconds. M ost theories of prospective tim- 
ing embody similar components!, including an internal time 
keeper, attention and memory23, A clock metaphor is used to 
describe the timekeeper mechanism, which represents subjective 
timethrough the accumulation or readout of pulses, possibly gen- 
erated by oscillators. Our perception of time, however, is inti- 
mately related to the level of attention given to the passage of time. 
When attention is diverted, a systematic shortening of subjective 
duration occurs, implying that pulses from the timekeeper may 
belost*. Attention may also mediate the flexible starting and stop- 
ping of pulses from the timekeeper, which enables anticipation of 
predictable events’. Hence, a representation of subjective time 
emerges from the interplay between timekeeping and attention 
mechanisms. This representation is then passed on to working 
memory, a short-term repository where interval representations 
are maintained and manipulated in accord with current goals (for 
example, comparing two intervals of time)® Working memory 
functions can therefore alter stored representations of time as well. 
The combination of these different component processes gives 
rise to the subjective perception of time, although the relative con- 
tribution of each might differ depending on the interval duration 
or the cognitive demands of timing events’. 

Theneural systems that support different component process- 
es of time perception are a matter of debate. The basal ganglia 
and lateral cerebellum have been logical candidates for hypo- 
thetical timekeeping operations, as damage to these brain regions 


nature neuroscience * volume 4 no 3 * march 2001 


commonly disrupts behaviors that depend upon precise timing, 
such as rhythmic movements in Parkinson's disease’ and regu- 
lation of agonist- antagonist muscle activity (for example, dys- 
metria) in cerebellar damage®. Although these movement 
abnormalities could be attributed to disruption of more gener- 
alized motor execution functions, the basal ganglia and cerebel- 
lum do seem to mediate time perception. Studies of Parkinson’s 
disease patients!" and pharmacological investigations in ani- 
mals!2:13 have argued that timekeeping operations are regulated 
through dopamine neurotransmission in the striatum. Human 
lesion studies indicate that the lateral cerebellar hemisphere and 
its primary output, the dentate nucleus!* 18 are also involved in 
timekeeping mechanisms. Nonetheless, it has been difficult to 
isolate timekeeping and attention operations from working- 
memory and response implementation processes!. Timing 
deficits after basal ganglia or cerebellar damage could also be due 
to abnormalities in interconnecting cortical systems commonly 
associated with some or all of these processes!9.20, Fewer studies 
have examined the involvement of the cerebral cortex in time 
perception. Focal lesion investigations in animals and humans 
have shown that the frontal and parietal lobes are also essential 
for accurate time perception, perhaps due to their purported 
attention and working memory functions!*21.22, Others have 
posited a role for the supplementary motor area2?, but this has 
been difficult to assess because focal lesions are uncommon in 
this region. 

Functional imaging techniques can be used to dissect the con- 
tribution of each component of multiple neural systems, 
although studies of timing using these methods have produced 
conflicting or ambiguous results to date’. M ost research2*-27 has 
focused on motor timing, making it difficult to separate activa- 
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Fig. 1. Trial events in the time perception (a), pitch perception (b), and 
control (c) conditions. In the time perception condition, subjects indi- 
cated whether the comparison interval (defined by tones 3 and 4) was 
longer or shorter than the standard interval (defined by tones 1 and 2). 
In the pitch perception condition, subjects indicated whether the com- 
parison tone (tone 4) was higher or lower in pitch than the standard 
tones (tones 1, 2 and 3). In the control condition, subjects pressed a 
key after the presentation of the four tones. 


tion in systems traditionally associated with motor control, such 
as the basal ganglia and cerebellum, from those supporting time- 
keeping or other cognitive processes. Two PET studies?829 have 
specifically examined time perception. Unfortunately, the time 
scale of PET scanning is limited to blocked-trial designs that can- 
not disentangle processing associated with encoding an interval 
from processing associated with decision making and imple- 
menting a response. We reasoned that fundamental insights into 
this issue could be gained by studying the time course of brain 
activation patterns associated with different components of a 
time perception task. The present study exploited the finer tem- 
poral resolution of event-related functional magnetic resonance 
imaging (fM RI) to isolate patterns of brain activation that cor- 
related with encoding time intervals from those associated with 
comparing two time intervals and implementing a response. Tim- 
ing theory suggests that activation in systems integrally involved 
in encoding or formulating arepresentation of time (pacemaker 
and attention operations) should develop at the onset of a to-be- 
timed event?3, followed by activation in systems concerned with 
manipulating information in working memory (comparing inter- 
vals) and implementing a response 

We obtained fM RI scans of seventeen subjects as they per- 
formed three different tasks, the order of which was counterbal- 
anced across subjects. In the time (T) discrimination condition, 
two tones (50 ms) separated by 1200 ms (standard tone pair) were 
presented, followed by a 1-s delay and then a comparison tone 
pair (Fig. 1a). Subjects indicated whether the comparison tone 
pair was longer or shorter than the standard. To better separate 
neural systems specific to timing, subjects also performed a pitch 
(P) discrimination condition in which the auditory events were 
similar except that subjects indicated whether the fourth tone was 
higher or lower in pitch than the first three tones (Fig. 1b). N eur- 
al systems involved with processing time and pitch information 
were identified by contrasting imaging runs in each discrimina- 
tion condition with a sensorimotor control (C) condition in which 
subjects responded after the presentation of two isochronous tone 
pairs of identical pitch (Fig. 1c). TheT and P conditions were then 
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contrasted to specify systems unique to time discriminations. 
These subtractions were conducted at each of four scanning inter- 
vals after trial onset (2.5, 5.0, 7.5 and 10.05). In all conditions, the 
typical motor response occurred approximately 4.5 s after trial 
onset (Fig. 2). Allowing 5 s for the hemodynamic response to peak, 
we proposed that the 2.5- and 5.0-s intervals after trial onset 
should reveal brain activation patterns specific to encoding time 
intervals. In contrast, the 10.0-s scanning interval should include 
activations associated with contrasting the standard and com- 
parison intervals and implementing the response. Overlap between 
these processes should be particularly evident during the 7.5-s 
scan, due to encoding of the comparison interval. The results 
reported here show early sustained activation of the basal ganglia 
and right inferior parietal cortex, implicating these systems in for- 
mulating representations of time. Though activation in the cere 
bellum was more robust during timethan pitch discriminations, 
activation was located in the vermis and unfolded late, suggest- 
ing amore general involvement in cognitive or sensorimotor func- 
tions. The evolution of activation in the bilateral prenotor and 
right DLPF cortex differed from each other, consistent with pre 
vious work implicating these systems in different aspects of work- 
ing memory. 
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Fig. 2. Temporal relationship among the trial events, acquisition of 
images and hypothetical hemodynamic response functions to different 
task components. Seven scans were acquired during each 17.5-s trial 
(a 2.5-s inter val between the seventh image and the first image of the 
next trial is not illustrated on the timeline). The first scan was 
acquired at the onset of the first tone (T1). The fourth tone (T4) was 
presented an average of 3.4 s after trial onset. The typical key press 
response occurred 4.5 s after trial onset. The two hypothetical time 
course functions illustrate early versus late MR signal responses to dif- 
ferent trial events. An early response corresponding with the encod- 
ing of temporal information (red plot) would have a maximal signal 
change at 2.5 and 5.0 s after trial onset. In contrast, a late response 
due to decision making and response preparation processes (blue 
plot) would be observed primarily at 7.5 and 10.0 s after trial onset. 
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RESULTS 
Behavioral data collected during scanning 
showed that response times and accuracy corre- 
lated with the difficulty of time and pitch dis- 
criminations. Reaction time was typically longer 
(Fig. 3a, F576 = 4.2, p< 0.01; Fig. 3c, F697 = 4.0, 
p < 0.01) and accuracy poorer (Fig. 3b, 
Fas7 = 8.1, p< 0.001; Fig. 3d, Fy 412 = 2.7, 
p < 0.025) when the comparison stimuli were 
closer in time or in pitch to the standard stimu- 
lus. There were no significant differences between 
the two discrimination conditions in overall 
accuracy (T, 83 + 3%; P, 78 + 3%) or reaction 
time (T, 1111 + 76 ms; P, 1076 + 54 ms). Reac- 
tion times for the C condition (707 + 39 ms) 
were significantly faster (F 11g = 48.9, p < 0.0001) 
than those for the time and pitch conditions. 
During the early imaging epochs (2.5 and 
5.0 s), which emphasize encoding of temporal 
information, subcortical activations specific to 
theT condition (Table 1) were observed within 
the right putamen, head of the caudate nucleus 
bilaterally, and right centromedian and ven- 
troanterior thalamic nuclei (Fig. 4a). Early acti- 
vation specific to the T condition was also 
observed in various cortical regions (Fig. 5): 
right intraparietal sulcus (BA 40), bilateral dor- 
sal and left ventral premotor areas (BA 6), and 
bilateral lateral temporal 
cortex (BA 21/22). Activa- 
tion specific to the T condi- 
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Fig. 3. Mean (+ standard error of mean) reaction time and percent correct for the time per- 
ception (a b) and the pitch perception (c, d) conditions. Data are depicted as a function of the 


comparison interval or comparison pitch. 


Table 1. Stereotaxic brain atlas coordinates”? for Time >C ontrol subtraction. 


tion was sustained during Location (Brodmann Area) Hemisphere 2.5 5.0 7.5 10.0 
the 7.5- and/or 10.0-s imag- Basal Gandli 
ing epochs in most of these Basal Ganglia 
regions. In theP condition, Medial caudate (head) h a : ` r cee 

r f activation durin Saera FA Basu 
a md -i n Eri caudate (body) Pp : 7 
overlapped with those in the utamen 1 Be 19 
T condition. In both the 26, 6, -2 
T and P conditions (Table 2), L -20, -1,5 
activity unfolded early with- Cerebellum 
in the medial wall (preSM A Vermis (tuber, lobule VIIB) B -3, -70, -30 2, -70, -29 
and SM A proper, BA 6 and Thalamus 
anterior cingulate, BA 32; Centromedian nucleus R 4, -21,0 4,-21,0 
Fig. 4c) and the anterior Ventroanterior nucleus R 4, -11,0 5,-10,0 
insula/frontal operculum Frontal 
(Fig. 4a), but was sustained Dorsal premotor (6) R 23, -7,48 23, -3, 52 46, 1, 49 
during later epochs as wall. L -45,-7, 47 

During the later imaging Ventral premotor (6) R 46,8, 24 
epochs (7.5 and 10.0 s), L -54,-13,26  -51,-15, 27 
which included decision and Dorsolateral (46/10/9) R 34, 23, 25 31, 46, 22 
response selection compo- 41,29, 22 
nents of the tasks, acti vation parietal 
specific to the T condition Int ideal 
(Table 1) was observed in eae 
the posterior vermis (tuber) Angular gyrus (40) R 38, -40, 41 36, -43, 40 37, -47, 38 30, -56, 35 
of lobule VIIB of the cere- Superior parietal lobule, 
bellum (Fig. 4b) and the Precuneus (7) R 10, - 68, 44 
right dorsolateral prefrontal © Temporal 
(DLPF) cortex (BA 46/10/9; Superior temporal (22) 51, -39, 6 
Middle temporal (21) L -46, -56, 4 


Fig. 5). All other activation 
foci were observed in the left 
hemisphere in both the 


R, right; L, left; B, bilateral. The activations reported in this table were not obser ved in the Pitch > Control subtraction. 
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T and P conditions (Table 2), and included the inferior frontal 
gyrus (Broca's area, BA 44/45), intraparietal sulcus (BA 40), supe 
rior parietal lobule/precuneus (BA 7) and DLPF cortex. 

The results from the T minus P subtraction were similar to 
the results for the T minus C subtraction (Fig. 6). During the 
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Fig. 4. Activation foci in the basal ganglia (a), cerebellum (b), and pre- 
supplementary motor area/anterior cingulate (c) resulting from subtrac- 
tion of the control (C) condition from the time (T) and the pitch (P) 
perception conditions at 2.5, 5.0, 7.5 and 10.0 s after trial onset. 
Significant foci (p < 0.001) are displayed with a red-yellow intensity scale 
denoting greater activation for the T or P conditions. Slices are displayed 
in neurological view (left is on the viewer's left). Location of slices defined 
by the distance (mm) from anterior commissure: x, right (+)/left (-); y, 
anterior (+)/posterior (-); z, superior (+)/inferior (-). Caud, caudate 
nucleus; Cing, anterior cingulate area; Ins, insula; O per, frontal opercu- 
lum; Put, putamen; Thal, thalamus; SMA, supplementary motor area. 


earlier imaging epochs (2.5 and 5.0 s), subcortical activations 
unique to the T condition werein the right hemisphere and 
included the putamen (x, y, z = 24, 7, -2), caudate (15, 6, 13) and 
insula/frontal operculum (29, 16, 2). The later region, however, 
was also activated during the 7.5-s epoch in the pitch condition 
(Table 2, Fig. 4a). During the later imaging epochs 
(7.55), theright DLPF cortex (21, 21, 30) was also uniqueto the 
T condition (Fig. 6). 


DISCUSSION 

The present findings provide compelling evidence for the involve- 
ment of the basal ganglia in formulating representations of time. 
Activation in the right putamen and caudate were uniquely asso- 
ciated with encoding time intervals. These results corroborate 
studies in Parkinson's disease showing that dopaminergic treat- 
ment improves motor timing?03! and time perception=2. Phar- 
macological challenges in animals also suggest that dopaminergic 
antagonists and agonists respectively slow down and speed up 
timing operations!?13, Contrary to one proposal?3, these and 
other studies!011.27 show that the basal ganglia are involved in 
timing a wide range of intervals, from hundreds of milliseconds 
(300 ms) to tens of seconds (20 s). Collectively, these results 
implicate striatal dopaminergic neurotransmission in hypothet- 
ical internal timekeeping mechanisms. 


Table 2. Stereotaxic brain atlas coordinates”? for regions commonly activated in subtractions of Time and Pitch 


perception conditions relative to C ontrol condition. 


Time >Control Pitch >Control 
Location (Brodmann Area) Hemisphere 2.5 5.0 10.0 2.5 5.0 75 10.0 
Frontal 
Insula/operculum (47) R 31,17,3 35, 16; 3 34,17,4 34, 17,0 
fi -35,11,5 -34, 15,2 -36, 12, 4 -34,18,1 -36, 17,0 
PreSMA (6), 
Anterior cingulate (32) L -4, -1,56 -4,6,49 -7,10, 45 -5, 12, 43 -6,7,48 -4,8, 49 
Inferior frontal gyrus (44/45) R 37,1, 32 37, 4, 28 
L -46, 4, 21 -47,5, 18 -45,4, 22 -44, 7, 26 
Dorsolateral (46/10/9) L -39, 42,12 -36, 46, 13 -36, 40,8 
-42, 26, 28 -40, 14, 29 
Parietal 
Intraparietal sulcus, 
Angular gyrus (40) L -31,-49,37  -29,-52, 33 
-36,-53, 44 -32, -47, 38-30, -55, 36 
Superior parietal lobule, 
Precuneus (7) L -21,-66,49 -28,-49, 43 -13,-72,50 -43,-57,50 
-21, -63,51 -25,-65, 50 


R, right; L, left; B, bilateral 
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Fig. 5. Activation foci in the lateral surface of the left and right hemi- 
spheres denote greater activation for the time (T) and the pitch (P) 
perception conditions relative to the control (C) condition at 2.5, 5.0, 
7.5 and 10.0 s after trial onset. Significant foci (p < 0.001) are displayed 
in red. DLPF, dorsal lateral prefrontal cortex; D. Premotor, dorsal pre- 
motor; IFG, inferior frontal gyrus; Ins, insula; IPS, inferior parietal sul- 
cus; MTG, middle temporal gyrus; Oper, frontal operculum; STG, 
superior temporal gyrus; V. Premotor, ventral premotor. 


Our findings did not support a unique rolefor the cerebdlum in 
encoding time intervals. Nonetheless, cerebellar activation was 
observed during the time perception task (T minus C), consistent 
with several studies showing diminished time perception in patients 
with cerebellar damage!®18.34, H owever, in our study, activation 
was in the vermis rather than the lateral cerebellar hemispheres, 
contrary to reports that damage to the lateral cerebellum, but not 
the vermis, correlated with time perception deficits!518, Cerebd- 
lar activation evolved later in the course of the trial, just before and 
during movement execution, suggesting an involvement in process- 
es other than explicit timing. This is consistent with our previous 
fM RI study?” showing that cerebellar activation was not specific to 
timing self-paced finger movements. Apart from its well-docu- 
mented rolein sensorimotor processing, neuroimaging research 
indicates that the cerebellum participates in many cognitive func- 
tions, including tactile perception? and working memory. One 
lesion study has also shown that cerebellar damage produces pitch 
perception deficits!4, Its broad role in sensorimotor and cognitive 
processing?’ has suggested that the cerebellum monitors and adjusts 
input from the cerebral cortex, but isnot involved in computing a 
specific operation per se*8, By this account, later activation in vermal 
lobule VIIB, which receives auditory and visual input32, could be 
due to its involvement in optimizing sensory input from auditory 
systems, which facilitates the comparison of intervals in working 
memory. Although other explanations are possible, this account is 
appealing because it predicts that damage to the cerebellum will 
slow sensory acquisition, which should disrupt a broad range of 
behaviors, especially those involving timing. This view may explain 
why patients with cerebellar damage show deficits in timing!©2’, 
but not always in the perception of pitch or loudness!®18, 

Representations of time depend on the interplay of internal 
timekeepers with attention and working memory, functions 
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Left hemisphere 


Right hemisphere 


more commonly identified with cortical systems. Neural sys- 
tems associated with these functions should support a variety of 
computations, which may explain why they were not always unique 
to timing intervals (T minus P). However, in the comparisons 
involving the control condition (T minus C, P minus C), right 
hemisphere activations were observed during time but not pitch 
perception. These later results are consistent with findings from 
converging neuroscience approaches. Specifically, a neu- 
roanatomical bridge for basal ganglia-cortical interactions is 
the thalamus”, which was activated early during the encoding 
of intervals, along with two cortical regions, suggesting they 
work together in formulating representations of time. Coupled 
activation in the right inferior parietal cortex may suggest an 
interdependent role of this region in attention, which theoret- 
ically regulates the timekeeping mechanism. Neurological 
patients with right but not left inferior parietal damage show 
time, but not pitch, perception deficits that correlate with 
impairments in switching attention2!. Electrophysio- 
logical recordings in humans have also shown aright 
hemisphere bias for temporal processing*!, especial- 
ly in the parietal cortex*2. The close relationship 
between timekeeping and attention is presumed by 
one influential theory2, and has received empirical 
support in behavioral studies conducted on 
humans*. According to this view, representations of 
time are reflected in the pulse count accumulated over 


Fig. 6. Activation foci in the basal ganglia, insula/frontal oper- 
culum and dorsal lateral prefrontal cortex resulting from 
greater activation for the time (T) relative to the pitch (P) 
perception conditions at 2.5, 5.0, 7.5 and 10.0 s after trial 
onset. Significant foci (p < 0.001) are displayed with ared-yel- 
low intensity scale. Slices are displayed in neurological view 
(left is on the viewer's left). Location of slices defined by the 
distance (mm) from anterior commissure. Caud, caudate 
nucleus; DLPF, dorsal lateral prefrontal cortex; Ins, insula; 
O per, frontal operculum; Put, putamen. 
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a particular physical time, which critically depends on the 
degree of attentional engagement. Our results point to the right 
inferior parietal cortex in regulating the accumulation of puls- 
es, because of its well-documented involvement in attention“. 
Bilateral projections from the inferior parietal cortex to the 
putamen and caudate nucleus in monkeys” provide a neu- 
roanatomic basis for the interaction of attention and time- 
keeping operations. 

The perception of time also relies on stored representations of 
intervals in working memory”. During time perception, activa- 
tion was observed in regions commonly associated with tempo- 
rary storage functions, including the bilateral premotor (BA 6) 
and right DLPF cortex (BA 9, 10, 46) 19:20.45, Right DLPF acti- 
vation was also uniqueto performing time discriminations. This 
corroborates our previous finding that damage to these same 
regions in the right, but not left, hemisphere produces time per- 
ception deficits?!. Controversy exists over whether these areas 
support different working memory functions*-47, However, a 
recent meta-analysis of neuroimaging studies% implicated the 
premotor cortex in a ‘rehearsal circuit’ in tasks involving main- 
ly the temporary maintenance of information, such as item 
recognition. In contrast, the DLPF cortex was associated with 
an ‘executive circuit’ in tasks requiring manipulation of stored 
information, such as the two- and three- back working-memory 
tasks. Our findings are compatible with this process distinction, 
as premotor cortex activation began early, consistent with the 
need for maintaining the standard interval during the trial, 
whereas DLPF cortex activation unfolded later in association 
with comparing the two intervals and selecting a response Inde 
pendent evidence for the DLPF cortex in executive functions of 
working memory was observed in the pitch condition as well, 
in which activation unfolded later during the comparison phase, 
but was confined to the left hemisphere. Though premotor cor- 
tex was not activated in the pitch condition, repeated presenta- 
tion of the standard pitch across the trial may have minimized 
the need for rehearsal. 

In summary, the present results are compatible with prevail- 
ing cognitive theory, and provide new insights into the evolu- 
tion of activation in cortical and subcortical systems that are 
specific to different cognitive components of a time perception 
task. The reciprocal interactions among these specialized sys- 
tems give riseto our perception of current time. The results are 
in agreement with converging avenues of research implicating 
a perceptual system in which the basal ganglia act as a timekeeper 
that is tightly coupled with an attention system in the right infe 
rior parietal cortex. This right hemisphere bias for the encod- 
ing of temporal information is in agreement with converging 
focal lesion and electrophysiological research in humans. The 
distinct evolution of activation in the bilateral premotor and 
right DLPF systems, together with previous neuroimaging stud- 
ies, provides evidence for different working memory functions 
underlying time perception. Our results also showed that time 
and pitch discriminations are mediated by shared parietal and 
prefrontal systems mostly in the left hemisphere, which were 
activated during decision and response selection components 
of both tasks. Presently, we are investigating the dynamics of 
brain activation patterns during longer delay periods to more 
directly distinguish systems involved in encoding and short-term 
maintenance of time intervals. 


METHODS 


Subjects. Right-handed subjects (2 male/15 female; mean age, 23.9 years) 
gave written informed consent and were compensated for participation. 
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The experimental protocol was approved by the institutional review board 
of the M edical College of Wisconsin. 


Experimental design. Tone stimuli were presented binaurally using a 
computer playback system. Sounds were amplified near the scanner and 
delivered to the subject via air conduction through 180-cm paired plas- 
tic tubes, which were threaded through tightly occlusive ear inserts that 
attenuated background scanner noise to approximately 75 dB sound pres- 
sure level (SPL). Background scanner noise consisted of pulses occur- 
ring every 205 ms throughout the imaging run; theintensity of the tone 
stimuli averaged 100 dB SPL. For all three conditions, the standard tones 
were 700 Hz in pitch separated by a 1200 ms interval 
(Fig. 1). In theT condition, the eight comparison intervals were +60-ms 
increments of the standard interval, and were presented twice in a ran- 
domized order (16 trials); pitch did not vary across the four tones. In 
the P condition, the eight comparison tone pitches were +4 Hz incre 
ments of the standard 700 Hz tones and were presented twice in a ran- 
domized order (16 trials); duration did not vary during this condition. 
In the C task, 16 trials of identical standard tones were presented. The 
C task was a baseline condition used for removing the effects of low- 
level sensory and motor processing from the functional maps in the two 
discrimination conditions. Subjects pressed one of two keys with thei 
right index or middle finger to indicate longer/higher or shorter/lower 
in the discrimination conditions; subjects pressed a key using their index 
in the control task. Accuracy and reaction time were measured with a 
nonferrous key- press pad. Subjects briefly practiced the three condi- 
tions before scanning. 


= 


Image acquisition. Event-related fM RI was doneon a1.5T GE Signa 
(Waukesha, Wisconsin) scanner equipped with a three axis local gradi- 
ent head coil and an elliptical endcapped quadrature radiofrequency coil. 
Foam padding limited head motion within the coil. Echo-planar images 
were collected using a single-shot, blipped gradient-echo echo-planar 
pulse sequence (TE, 40 ms; TR, 2.5 s; 90° flip angle FOV, 240 mm; reso- 
lution, 64 x 64 matrix). Seventeen contiguous sagittal 7-mm-thick slices 
were acquired to provide coverage of the entire brain. Scanning was syn- 
chronized with the onset of the first tone so that 7 images were acquired 
during each 17.5-s trial (Fig. 2) with a total of 112 images per run (16 tri- 
als per run). An additional 4 images (10.0 s) were added to the beginning 
of the run to allow the M R signal to reach equilibrium, and were discarded 
from further analysis; 4 images were added to the end of the run to accom- 
modate the delayed rise of the hemodynamic response. Before function- 
al imaging, high-resolution three-dimensional spoiled gradient-recalled at 
steady- state anatomic images were collected (TE, 5 ms; TR, 24 ms; 40° 
flip angle; N EX, 1; slice thickness, 1.2 mm; FOV, 24 cm; resolution, 256 
x 128) for anatomic localization and co-registration. 


fMRI data analysis. Functional images were generated using Analysis of 
Functional N eurol mages“8 software. Time series images were spatially reg- 
istered in three-dimensional space to minimize effects of head motion. A 
deconvolution analysis was used to generate impulse response functions 
(IRFs) of thefM RI signal on a voxel-wise basis. This analysis produced an 
estimate of the hemodynamic response for each condition (T, P and C) rel- 
ative to a baseline state (rest) without making a priori assumptions regard- 
ing the shape, delay or magnitude of the IRF. Anatomical and functional 
images were then interpolated to volumes with 1 mm: voxels, co-registered, 
converted to Talairach stereotaxic coordinate space’, and blurred using a 
4mm Gaussian full-width half-maximum filter. Voxel- wise analyses of vari- 
ance (T versus C, P versus C, and T versus P conditions) were done sepa- 
rately for images obtained at 2.5, 5.0, 7.5 and 10.0 s after trial onset. 
Pooled- variance t-tests were applied on a voxel-wise basis to the IRF esti- 
mates for each epoch to identify regions showing greater activation in theT 
and P discrimination conditions relative to the C condition and greater acti- 
vation in theT than the P condition. An activated region was defined by 
an individual voxel probability less than 0.001 (t > 3.61; df, 16), and a min- 
imum cluster size threshold of 300 microliters. These two thresholds were 
established based on 10,000 M onteCarlo simulations demonstrating that the 
chance probability of obtaining a significant activation cluster for an entire 
volume (typel error) was less than 10-6. 
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